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ABSTRACT 

We compute both extrinsic (lensing) and intrinsic contributions to the (galaxy- )density-ellipticity cor- 
relation function, the latter done using current analytic theories of tidal alignment. The gravitational 
lensing contribution has two components: one is analogous to galaxy-galaxy lensing and the other arises 
from magnification bias - that gravitational lensing induces a modulation of the galaxy density as well 
as ellipticity. On the other hand, the intrinsic alignment contribution vanishes, even after taking into 
account source clustering corrections, which suggests the density-ellipticity correlation might be an in- 
teresting diagnostic in differentiating between intrinsic and extrinsic alignments. However, an important 
assumption, commonly adopted by current analytic alignment theories, is the Gaussianity of the tidal 
field. Inevitable non-Gaussian fluctuations from gravitational instability induces a non-zero intrinsic 
density-ellipticity correlation, which we estimate. We also argue that non-Gaussian contributions to the 
intrinsic ellipticity- ellipticity correlation are often non-negligible. This leads to a linear rather than, as 
is commonly assumed, quadratic scaling with the power spectrum on sufficiently large scales. Finally, 
we estimate the contribution of intrinsic alignment to low redshift galaxy-galaxy lensing measurements 
(e.g. SDSS), due to the partial overlap between foreground and background galaxies: the intrinsic 
contamination is about 10 — 30% at 10'. Uncertainties in this estimate are discussed. 
Subject headings: cosmology: theory — gravitational lensing — large-scale structure of universe — 
galaxies: halos — galaxies: structure 



1. INTRODUCTION 

Weak gravitational lensing provides a direct probe of 
the large scale distribution of mass (Gunn 1967,Miralda- 
Escude 1991, Blandford et al. 1991, Kaiser 1992; see Mel- 
lier 1999 and Bartelmann & Schneider 2001 for reviews). 
Analyses by several groups have demonstrated that mea- 
surements of galaxy shapes can be made sufficiently pre- 
cise to detect such an effect (Wittman et al. 2000, van 
Waerbeke et al. 2000, Bacon et al. 2000, Kaiser et al. 

2000, Maoli et al. 2001, Rhodes et al. 2001). As the 
precision improves, it is important to understand possible 
small contaminations to the lensing signal. An assump- 
tion that has generally been made is that there is no in- 
trinsic correlation between shapes/orientations of separate 
galaxies, other than through gravitational lensing (extrin- 
sic alignment). Recently, it has been pointed out that 
some amount of intrinsic alignment might be expected due 
to long range tidal correlations (Lee & Pen 2000, Pen, Lee 
& Seljak 2000, Croft & Metzler 2000, Heavens et al. 2000, 
Brown et al. 2000, Catelan et al. 2001 Crittenden et al. 

2001, Mackey et al. 2001, Porciani et al. 2001, Vitvitska 
et al. 2001, Mailer et al. 2001, van den Bosch et al. 2002 
) . The estimated level of contamination ranges from weak 
(< 10%) for source galaxies at a median redshift of 1, to 
dominant for galaxies at a redshift of 0.1 or smaller. 

There are several ways one might be able to differenti- 
ate intrinsic from extrinsic correlations. This is important 
not only from the point of view of isolating the lensing 
signal; the intrinsic correlation signal is also interesting in 
its own right - unraveling its origin may teach us much 
about the origin of angular momentum of galaxies. First 
of all, if a sufficient number of galaxies with accurate red- 
shifts (say using photometric redshifts) is present, one can 



investigate the scaling of the measured correlations with 
redshift - the intrinsic and extrinsic signals are expected 
to vary differently with the median redshift of the source 
galaxies, as well as the width of the galaxy distribution. 
At the moment, such an approach might be challenging 
to implement due to the lack of a large number of accu- 
rate redshifts in typical deep lensing surveys. Second, as is 
well-known, gravitational lensing introduces an ellipticity- 
ellipticity correlation that is curl-free (e.g. Stebbins 1996) 
n , while as emphasized by Pen et al. (2000) and Crit- 
tenden et al. (2000), intrinsic correlations generally carry 
both a curl and a curl-free component (or a magnetic and 
an electric part). Pen, van Waerbeke & Mellier (2001) 
(see also Hoekstra, Yee & Gladders 2002) presented some 
recent measurements. However, the curl part of the mea- 
sured ellipticity correlations is also often used as a diagnos- 
tic for systematic errors such as those introduced by corre- 
lated point-spread- function variations. It would therefore 
be useful to study other statistics that might help separate 
the intrinsic from the extrinsic signals. 

This was what motivated us to study the density- 
ellipticity cross correlation. Most studies so far have 
focused on ellipticity-ellipticity correlations (ei(9)ej(9')) , 
where i and j stands for the two components of the ellip- 
ticity, and 6 and 6' denote two angular positions in the 
sky. What we would like to investigate in this paper is the 
density-ellipticity correlation (S g (9)ei(6')) , where 5 g (0) is 
the galaxy overdensity (<5 S = (n g — n g )/n g where n g is the 
galaxy number density, and n g is its mean) at angular po- 

1 Schneider et al. (2001) recently pointed out there is a non- 
zero curl contribution (B-mode) from lensing due to source redshift 
clustering. The signal is only significant, however, at small angular 
scales, below about an arcminute or so. 



1 



2 



sition 0. ^ This quantity is very much analogous to what 
is often measured in the context of galaxy-galaxy lensing 
(Brainerd et al. 1996, Schneider 1998, Fischer et al. 2000). 
The difference is that while in galaxy-galaxy lensing d g is 
measured from a foreground population of galaxies and 
from a background population, here, we are mainly in- 
terested in 5 g and £j measured from the same population 
of galaxies (although we will have the occasion to study 
galaxy-galaxy lensing as well). 

What makes the density-ellipticity correlation, (S g ei), 
particularly interesting is the fact that intrinsic alignment 
makes no contributions to it, at least according to cur- 
rent tidal theories of intrinsic alignment (Crittenden et al. 
2001 [CNPT hereafter] Mackey et al. 2001 [MWK here- 
after]). This is easy to see given that S g oc V 2 </> (where <fi is 
the gravitational potential, and linear biasing is assumed, 
which should hold on large scales, where current tidal the- 
ories are supposed to work), and the intrinsic ellipticity 
depends quadratically on <fr. Gaussianity of the gravita- 
tional potential fluctuation, assumed by both CNPT and 
MWK, guarantees that the expectation value of the re- 
sulting cubic product of <f> from (S g ei) vanishes. What is 
less obvious is whether this continues to hold if the source 
clustering correction is included (which we label as (S c 
and scales quartically with < 
that it does. 



We will demonstrate in § 2.1 



In §2.2, we work out the gravitational lensing (extrin- 
sic) contribution to the density-ellipticity correlation. It 
turns out to have two components, one very much analo- 
gous to galaxy-galaxy lensing, and the other arising from 
magnification bias. The relative importance of the two de- 
pends on the redshift distribution of the galaxies, which 
we illustrate with a few examples. 

It appears then that the density-ellipticity correlation 
might provide an interesting test of the lensing hypoth- 
esis in ellipticity- ellipt icity measurements, particularly in 
view of the fact that intrinsic alignment creates no density- 
ellipticity correlation according to current tidal theories. 
However, we expect the latter to break down when one 
takes into account inevitable non-Gaussian fluctuations in- 
troduced by gravitational instability. This is discussed in 
§||. In fact, we will argue that the existing calculations of 
the intrinsic ellipticity- ellipticity correlation, by ignoring 
non-Gaussian contributions, miss a term that scales lin- 
early with the power spectrum, and that therefore domi- 
nates over the Gaussian (quadratic) term on large scales. 
A simple but crude way to understa nd the origin of this 
result is e ncapsulated in eq. 
In 
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we apply the non-Gaussian calculation to esti- 
mate the contamination from intrinsic alignment in galaxy- 
galaxy lensing measurements. As we have discussed, 
galaxy-galaxy lensing measures the same quantity (S g ei) 
except that S g is computed from a foreground population 
and €i from a background population of galaxies. Because 
this separation into two populations is usually imperfect, 
the overla p all ows some contamination from intrinsic align- 
ment. In §3.4, we discuss briefly the inclusion of a nonlin- 
ear galaxy-bias. Finally, we summarize our findings in §|| 
A discussion of the estimators for the density-ellipticity 
correlation, in both real and Fourier space, can be found 

2 We will also examine a correction to this quantity due to source 
clustering: <5 g (0)ef (0')); its precise definition will be given in §bj. 



in the Appendix. 

A reading suggestion: for readers interested in the non- 
Gaussian contributions to the intrinsic cllipticity-cllipticity 
and density-ellipticity correlations, they can skip directly 
to §[| Much of §|| describes a derivation that, under the 
assumption of Gaussianity and linear bias, current theories 
imply a vanishing intrinsic density-ellipticity correlation - 
parts of this are somewhat complicated, e spec ially those 
concerning source clustering corrections (§2.1), and can 
be skipped by readers not interested in the details; the 
beginning of §[^ is still recommended for introducing our 
notations. 

Before we begin developing our formalism and defining 
our notations in §^, it is helpful to mention some pre- 
vious work on related subjects. Kaiser (1992) discussed 
the density-ellipticity correlation in the context of weak- 
lensing. His calculation ignored the contribution from 
magnification bias. Lee & Pen (2001) discussed galaxy 
spin-density correlation, but their focus was on a quantity 
different from the ones we consider here, and their empha- 
sis is on application to a three-dimensional survey. That 
density-ellipticity correlation might be a useful quantity 
to consider in the context of both intrinsic and extrinsic 
alignments was mentioned by Catelan et al. (2001) but 
no explicit calculations were carried out. Catelan & Por- 
ciani (2001) discussed the density-tidal field correlations, 
but not the density-ellipticity correlation. 

2. DENSITY-ELLIPTICITY CORRELATION 

The ellipticity of a galaxy can be defined using the 
quadrupole moment Qij of the light distribution: Qij = 
J d 2 99i9jl(0), where the origin is chosen to be the centroid 
of the image, and 1(0) is the intensity profile. 
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This definition of ellipticity is also consistent with e\ = 
e cos 2a and €2 = e sin 2a, where e = (a 2 — b 2 )/(a 2 + b 2 ), 
with a and b being the major and minor axes and a being 
the orientation angle of the major axis. 

It is customary to consider the following two different 
ellipticity-ellipticity correlations (or their Fourier trans- 
forms, or the corresponding variance). Suppose A9 — A8 
(cos ipe, smtpg) is the angular separation vector between 
two galaxies (or two pixels where ellipticity estimates ex- 
ist), with ipg being the orientation angle between the sepa- 
ration vector and one's chosen x-axis, then one can define 
the tangential and radial ellipticities: 



e t = ei cos 2ipg + e 2 sin 2tp e 
e r = —ei sin 2%j)g + £2 cos 2ipg 



(2) 



which are equivalent to the ellipticities (ei and e 2 ) if 
the x-axis were chosen to lie along the separation vec- 
tor. The two corresponding correlations are (e t (6)e t (6')) 
and (e r (6)e r (6')), which depend only on the separation 
A8 = \9 — 9'\ but not its orientation. The cross-correlation 
between tangential and radial ellipticities vanishes by par- 
ity invariance. The electric and magnetic correlation func- 
tions can be constructed from combinations of the above 
correlations and their derivatives (see Kamionkowski, 
Kosowsky & Stebbins 1997, Seljak & Zaldarriaga 1997). 
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We are interested, on the other hand, in the density- 
cllipticity correlation: (S g (8)e t (8')), where S g is the galaxy 
overdensity. The correlation between density and the ra- 
dial ellipticity vanishes by parity. 

The observed ellipticity is divided into intrinsic and ex- 
trinsic parts: 

ti = + (3) 

where i can stand for 1 and 2, or t and r. 

Note that all of the above quantities are observed in 
projection i.e. 

e0) = J d Xg W g (x g )ei(e,x g ) (4) 

where Xg is the comoving distance along the line of sight, 
and we use e; with an extra argument of x§ to denote 
the ellipticity at a given angular position and a particular 
redshift (and similarly for ef 1 - and e° x ). Here W g (x g ) rep- 
resents the distribution of galaxies along the line of sight, 
whose normalization is J dx g W g (x g ) = 1. We will refer to 
it as the selection function. 

The galaxy overdensity is similarly projected [J 

5 g (d) = J dXgW g (Xg)S g (9,Xg) (5) 

The metric we adopt is: 

ds 2 = -dt 2 + a(t) 2 [d X 2 + r(x) 2 (d9 2 + sm 2 6d<p 2 )} (6) 

where r( X ) = K' 1 / 2 sinK^x, {-K)- 1 / 2 sinh(-/sT) 1 / 2 x, 
X for a closed, open and flat universe respectively, and 
K = —flkH^/c 2 with H being the Hubble constant to- 
day, c is the speed of light and Ofc is the fraction of critical 
density in curvature; a(t) is the Hubble scale factor as a 
function of proper time t. In this paper, we will display ex- 
amples exclusively for the model where the matter density 
Q m = 0.3, and the vacuum density Q\ = 0.7, although 
the formalism allows other possibilities (see e.g. Hui 1999, 
Benabed & Bernardeau 2001, Cooray & Hu 2001, Huterer 
2001, Hu 2002 for studies of quintessence models in the 
lensing context). 

Before we start describing our calculations in detail, it is 
important to emphasize that because the ellipticity et can 
only be estimated where there are galaxies, the ellipticity 
is always implicitly weighed by the local galaxy density. In 
other words, the actual quantity that can be realistically 
measured is: 

£cross(|0 - = J dXgW g ( X g) J d X ' g W g ( X ' g ) (7) 

(6 g (o,x g )(i + 6 g (e',x' g ))et(e',x! g )) 

For convenience, we define the following quantity: 

e t A (0) = J dXgW g ( X g)5 g (8,Xg)e t (8,Xg) (8) 

Then, Ccross = (5 g (8)e t (8')} + (5 g (8)ef(8')), instead of 
simply (5 g (8)et(0 )}■ This sort of correction is sometimes 
referred to as a source clustering effect (see Bernardeau 

3 This expression should in fact be_alightly modified by magnifi- 
cation bias, which we will discuss in §|2,2|, 



1998 who discussed this effect in the context of extrinsic 
ellipticity-ellipticity correlations). To be concrete, an esti- 
mator for ^ cr oss is given in the Appendix. We also provide 
an estimator for its Fourier analog. 

In this paper, we will ignore the source clustering cor- 
rection in the context of lensing (as is commonly done for 
ellipticity-ellipticity correlations). On the other hand, as 
far as intrinsic alignment is concerned, since, as we will 
see, the lowest order term {8 g e t ) vanishes (according to 
Gaussian theories), we will consider the term (S g e^) as 
well. We will use the symbol e Aln ' to denote the intrinsic 
part of 

2.1. Intrinsic Alignment - Gaussian Theories 

We calculate here both (Sgcf 1 ') and (<5 g e Ain ') using two 
different formulations of tidal alignment theories. 

2.1.1. Tidal Alignment Theory According to MWK 

We start with the simpler theory, developed by MWK 
(see also Catelan et al. 2001), who postulated that 

4 n - = 4[(2033 ~ 011 ~ 02 2 )(011 ~ 02 2 ) (9) 

15 

+3(^3 - 013)] 

e f' = T7 [2(2033 - 011 - 022)012 - 60i30 23 ] 

15 

where fj is a constant which quantifies the mean ellipticity 
of galaxies. The above is motivated by the tidal torquing 
theory for the origin of angular momentum developed by 
Peebles (1969), Doroshkevich (1970) and White (1984), 
and assuming the moment of inertia tensor is uncorrelated 
with the local tidal field. Here, the third direction is taken 
to be along the line of sight, is the gravitational potential 
at some given point in space, and 0y- = ViVj0. 

We relate 5 g to using the linear bias model and the 
Poisson equation i.e. 

S g = bS (10) 
V 2 = 4>rrGpa 2 S = (3Hfo m /2/a)8 

where 8 is the mass overdensity (5 = (p— p)/p with p being 
the mass density and p its mean), and b is the bias factor, 
which can be redshift dependent in general (Kaiser 1984, 
Bardeen et al. 1986, Fry & Gaztanaga 1993, Mo & White 
1996). g 

Limber's approximation, which assumes that the selec- 
tion function is slowly varying compared to the correlation 
function (Peebles 1980), tells us how to relate projected 
correlations to their three-dimensional counterparts: 

(* fl (0)e<(0')) - J d Xg W 2 J dpsiSgfateiq)) (11) 

(6 g (0)e? (6')) = J d Xg W 2 J dp 3 (5g(p)e i (q)5 g (q)) 

4 Recent more sophisticated modeling of galaxy biasing using the 
halo model still respects linearity on large scales (see e.g. Peacock & 
Smith 2000, Seljak 2000, Ma & Fry 2000, Scoccimarro et al. 2001, 
Bcrlind & Weinberg 2001). See also Tegmark et al. (2001) for a 
measurement of the galaxy-mass correlation coefficient, which is con- 
sistent with a deterministic linear bias on large scales. Note that the 
tidal theories proposed by MWK and CNPT are also expected to be 
valid only on large scales. 
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where the selection function W g (Xg) is as defined in eq. 
(H) & (||), P and q arc three-dimensional positions, the 
z-axis (labeled 3) is the direction along the line of sight, 
and p± — q^ (components of p and q perpendicular to z) 
is chosen such that pj_ — q_i_ = A8r(x g )- 
From eq. ( |To| ) and (^|), it is easy to see that 

(<Kp)ef (q)) = £ c AS (V 2 0(p)0 A (q)0 B (q)) (12) 

A,B 

where A and B signify various second derivatives, and cab 
represents coefficients independent of spatial positions (ex- 
cept through the overall redshift). That has equal prob- 
ability of being positive or negative implies the above must 
vanish (for instance, for a Gaussian random distribution 
of 0, as is assumed by MWK; see §|| for a discussion of 
expected violations of this assumption). That the above 
vanishes for both i = 1 and i = 2 implies (6 g (0)et(9')) 
must be zero as well when a projection is applied. 

How about (5 g (p)e- n (q)5 s (q)), which shows up in the 
source clustering correction (eq. JIM)? This requires a 
little bit more work. From eq. (Bjjit is clear there are 
several terms. We will work out the case of i = 2 in detail. 
The case of i = 1 is very similar, albeit with more terms. 



We have (6 g e 
given by 



<;'2 Sg) 



= b (Se™'5), and the latter factor is 



(5(p)e 2 n -(q)%)) 



2/3 
15 



[2<<5(p)<5(q)<Mq)0i2(q)) (13) 

-<s(p)<Kq)<Mq)<Mq)) 

-(<5(p)<5(q)0 22 (q)0i 2 (q)) 
-3(S(p)%)0i 3 (q)023(q)>] 



By MWK's assumption of a Gaussian random field 0, 
each of the above terms factorizes into products of second 
moments. For instance, the first term gives 

2(5(p)«5(q)033(q)0i 2 (q)) =2(5(p)033(q))(5(q)0i2(q)) (14) 

+2<5(p)<Mq)X%)<Mq)) 
+2<<5(p)<5(q)><<Mq)<Mq)) 

Several of the terms above contain an average of the 
following form: 



Wp)<Mq)> 



i 



/ -0f p 4k)k 2 hk je -*^ (15) 



inGpa 2 J (2tt) 



where is the potential power spectrum in the sense that 
(0(p)0(q)) = J(d 3 fc/[27r] 3 )P (fe) exp [-ik • (p - q)]. 

The above expression can also be used to evaluate 
(<5(q)0y (q)) by setting p = q. In fact, doing so, it is 
clear that (#(q)0i2(q)) vanishes by isotropy (i.e. P<j>(k) 
depends on magnitude of k but not its direction). We can 
therefore ignore the first term on the right hand side of eq. 

©• . 

Similarly, a term like (033(q)0i2(q)) must vanish, be- 
cause this is given by 



(033(q)^ia(q)> = 



d 3 k 



(16) 



Therefore, the only term that would eventually survive 
from eq. © is 2(<5(p)0 12 (q))(<5(q)0 3 3(q)}. This is (ig- 
noring the factor of 2/3/15) from the first term on the 
right side of eq. (p"3|). Most of the rest of the terms 
from eq. (|l^) vanish due to the same reasons as above 
except for two terms (which come from the second and 
third terms in eq. M): -(£(p)0i2(q))(£(q)0ii(q)) - 
(S (p) 0i2 (q) ) {5 (q) 022 (q)) ■ Finally, collecting all the sur- 
viving terms, it is easy to see that 



2(«5(p)0 12 (q))(5(q)033(q)) 
-<5(p)0i 2 (q))(5(q)0n(q)) 

-(<s(p)<Mq)>(<s(q)& 2 (q)) = o 



(17) 



because (5(q)03 3 (q)) = (<5(q)0 22 (q)) = (<5(q)0n(q)). 
This follows again from isotropy. 

Therefore, the conclusion is that (£(p)e 2 n - (q)<5(q)) as 
given in eq. ( Jl3| ) vanishes exactly. The same statement 
can be proven Tor (<5(p)e I 1 n (q)i5(q)), using very similar ar- 
guments as those outlined above. Hence, in summary, we 
find that (using the above together with eq. [pT|) accord- 
ing to the tidal theory of MWK, and under the assumption 
of a linear bias, both (6 g (0)ei(0')) and its source clustering 
correction (6 g (8)e^(8')) vanish exactly for both i = 1 and 
i = 2, and therefore, for i — t as well. 

2.1.2. Tidal Alignment Theory According to CNPT 

Next, we turn our attention to the tidal theory of CNPT 
(which built upon earlier work by Lee & Pen 2000), accord- 
ing to which the galaxy ellipticity is given by: 



2 ^2(4>lihl - <fe<fe) 
i=l 
3 

1=1 



(18) 



which vanishes by isotropy. 



where a is a constant which depends on both the average 
ellipticity of galaxies and the degree to which the moment 
of inertia tensor is correlated with the local tidal field 0, 
(f>ij is the unit normalized traceless tidal tensor i.e. 0y = 
[ViVj0 — ^5ijV 2 <f)]/N, where N is a normalizing factor 
so that j i'ij'fiij — 1 (summation of i, j is over 1,2,3). 
As such, this theory relates the galaxy ellipticity to the 
direction of the galaxy angular momentum, but not its 
amplitude. |^] 

Here, we would like to evaluate the density-ellipticity 
correlation using Limber's approximation and linear bias 
just as we have done using the other tidal theory (eq. 
Q & p|). We therefore, need to compute (5(p)ef'(q)) 
and ((5(p)e- n (q)(5(q)) at two spatial points p and q just 
as before. To do so, we need the two-point probabil- 
ity distribution for 0y. Adopting the notation T = 
(0n,022,033,0i2,0i3,023) where T signifies the whole 
vector of different components of the tidal field (i.e. T\ = 

5 The symbol a here is actually equal to aa(6 — 9iz/4) in the 
notation of CNPT. 

e Porciani et al. (2001a, b) recently showed using numerical sim- 
ulations that much of the angular momentum correlation in dark 
matter halos induced by tides is erased by non-linear effects. One 
can think of modeling this by decreasing a in eq. (tq). See also van 
den Bosch et al. (2002). 



5 



0n, etc.), the two-point probability distribution used by 
CNPT is Gaussian random i.e. 



P(T(p),T(q)) 



1 



\/dctC(27r) 6 



exp 



1 



fT(j-lf 



where P(T(p), T(q))d 6 T(p)d 6 T(q) gives the probability 
that at spatial points p and q, the tidal field vectors take 
their respective values in the above ranges. Here T = 
(T(p),T(q)), and T T is its transpose. The matrix C gives 
the correlation matrix, which has the block diagonal form: 



C 



Co C, 



Cp-q 



p-q 



C 



(19) 



where Co is the 6x6 zero-lag correlation matrix, and C p _ q 
is the 6x6 two-point correlation matrix. 
Let us first evaluate (<5ef ')■ 



■(q)) = E^ B (V 2 0(p)<Mq)0 B (q)) (20) 
A.B 



where A and B denote appropriate double indices as given 
in eq. (|l8|), and cIab represents coefficients independent 
of spatial positions (except through the overall redshift). 

Clearly, the probability distribution in eq. ( p^ ) is invari- 
ant under T — ► —T (i.e. <j) — ► — 0), while the combination 

V 2 4>4>a4>b (eq. [p0[ ) switches sign under such a trans- 
formation. Therefore, the expectation value (<5(p)eJ n - (q)) 
must vanish. 

Next, let us consider: 

(5(p)ef-(q)5( q )) (21) 

= J2 eA B (V 2 </)(p)^(q)0 B (q)V 2 0(q)) 

= ^e AB (TrT(p)TrT( q )f j4 (q)f s (q)) 

A,B 

where f = (0u, 022, 033, 012, 013, 02s), TrT = Ti + 
T2 + T3, and cab represents some coefficient. Follow- 
ing CNPT, it is useful to rotate the vector T by defining 
T = RT, where R is an invertible matrix such that T = 
(TrT/V3, (Ti - T 2 )/V2, (Ti +T 2 - 2T 3 )/v / 6, V2T 4 , y/2T 5 , 
V2Tq). With this rotation, the above expression can be 
rewritten as 



WpK"-(q)<5(q)} 

= £ e ^ (?i (P)3i (q)^4 (q)f fl (q)) 



(22) 



A.B 



where e' AB — 

Ed E e DER D 1 A RE B , and f = RT. Note 
that because T is formed from the traceless part of the 
tidal tensor, Ti = be definition. It is also useful to 
define |T| = [T 2 2 + T 3 2 + ... + T 6 2 ] 1 /2 j so that T A = f A \T\ 
for A^l. 

The advantage of such a rotation is that the correlation 
matrix at zero-lag becomes diagonalized: Cq = RCqR t is 
diagonal: Co = (£o/15) diag (5, 2, 2, 2, 2, 2), where £0 is the 
zero-lag correlation of TrT. Similarly, we define C p _ q = 



RC p -qR T , and 



Co 



°p-q 
Co 



(23) 



but note that neither C p _ q nor C is diagonal in general. 

The expectation value {T\TiTaTb) can be written as 
follows. 



(Ti(p)Ti(q)t 1 (q)t B (q)) 
d 6 T(p)d 6 T(q) 



(24) 



\/detC(27r) 6 
I 



T 1 (p)r 1 (q)T A (q)T B (q) 



x exp 



--T T C- l T 
2 



d 6 r( P )d 6 r(q) 

%/dctC(27r) 6 



Xi(q)T B (q) 



9 



SCfy 1 



— rexp 



-if^T 

2 



<9C 



T <T A (q)T B (q)) 
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a(detC)- 1 / 2 



V^detC(f A (q)T B (q)) 



where C-^ 1 represents the (1, 7) component of the symmet- 
ric matrix C. Note that d 6 T = dTr T\T\ 4 d\T\d 4 f /y/3. 

The second term in the last equality above (the one in- 
volving a derivative of the determinant), when contracted 
with e' A B m ec i- (Hl)i gives us something that is propor- 
tional to (e™'(q)), which must vanish by isotropy. There- 
fore, all we have to worry about is the first term. At 
first sight, it appears that the same argument might ap- 
ply: after contracting it with e' AB , the first term gives a 
derivative of (e- n -) with respect to C^ 1 . Should this also 
vanish by isotropy? The ultimate answer will turn out to 
be yes, but one has to reason with some care: does the 
fact that we are varying C^ 1 while keeping the rest of the 
components of C fixed somehow spoil isotropy? One would 
expect not, because the component we are varying has to 
do with the trace of the tidal tensor. However, let us prove 
this explicitly. 

The matrix C _1 in general takes the following form: 



c- x = 



M 



M 



where M is some matrix that satisfies 



C C~^C M - 



C p _ q M 



(25) 



(26) 



to ensure CC _1 gives the identity. 

When we are varying C-j^. 1 , we are therefore varying 
M11 while keeping everything else fixed. This means that 
A(Cpi q C Af) = 0. Eq. © should remain valid under 

the variation, and so -A(C Cpi q C M) + A(C p _ q M) = 0. 
Combining both, together with eq. (p6|), it can be shown 
that varying Ciy 1 implies a change in Co of the form 



AC = -C p _ q AMCo - CoAMCp- c 



(27) 



The fact that the only non-zero component of AM is 
A M11 = AC-,^. 1 , together with the fact that Co is diago- 
nal, implies that [ACo],j = except for i or j (or both) 
= 1. This is an important fact we will make use of later. 
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Dropping the q label for simplicity, let us then consider 
the term from eq. 



d 



oc 



— {TaT b 



(28) 
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d f d 6 T ~~ 

r / — . TaJb exp 

dCrf J \/dcTC^(2^) 3 

r d G T 



VdetC (27r) 



T A T B 



\t t C^AC q C^T 



exp 



where Q is a term that involves the derivative of the deter- 
minant det Co with respect to C^j 1 , and Q is proportional 

to (TaTb), which as we have argued before, vanishes by 
isotropy under contraction with e' AB (eq. p2| ) because it 
gives something that is proportional to (e 1 / 1 ). 

Now, recall the fact that (ACo)ij is non-zero only if i 
or j = 1. This, together with the fact that Cq 1 is di- 
agonal, implies that aside from Q (which we can ignore), 
all terms in eq. ( pgj ) are of the form t](TaTb{Ti) 2 ) or 
J2c Cci^A T'bT\Tc\T\) where r\ and (c are some coeffi- 
cients. The fact that the latter term involves an odd num- 
ber of directional vectors T implies it must be zero. 

Therefore, the only term we need to consider is 

(<Kp)er(q)%)) =^^ s (f A (q)f s (q)[Trr(q)] 2 ) (29) 

A,B 

= J2veAB(f A (ci)fB(ci)[Ti'T( q )} 2 ) 

A,B 

Note that r\ is in principle dependent on p — q. The co- 
efficients cab are determined by the expressions for ef 1 ' 
in eq. (|l8|). For i = 1, the above is proportional to 
((T ifi - T 2 T 2 ) ( Tr T) 2 ) + {(Tufa - f 6 f 6 ) (TrT) 2 ).^ The im- 
portant point to keep in mind is that this expectation value 
is evaluated at a single point (q), and with no preferred 
direction, this clearly vanishes by isotropy. Similarly, the 
i = 2 term vanishes as well. 

To summarize, we find that according to the tidal theory 
formulated by CNPT, and under the assumption of linear 
bias, both {8 g (B)ei{9')) and its source clustering correction 
(S g (8)e^(8')) vanish exactly. This holds for i = 1, i = 2 
and therefore for i = t as well. 

2.2. Extrinsic Alignment - Gravitational Lensing 

Gravitational lensing induces a correlation between el- 
lipticity and galaxy density, which we will work out in this 
section. The projected ellipticity and galaxy density are 



given in eq. (|4|) and eq. (I 
to be slightly modified to ta 



). The latter equation needs 
le into account magnification 



bias (e.g. Broadhurst, Taylor & Peacock 1995; Moessner, 
Jain & Villumsen 1998): 

6 g (0) = J dx g W g (x g )[5 a (e, Xg ) + 5(s-0A) K (e, Xg )} (30) 

where s is the slope of the luminosity function in the fol- 
lowing sense: if N(m, z) gives the surface number den- 
sity of galaxies per unit magnitude m per unit redshift z, 



dlogN(m, z)jdm = s. Note that s can be redshift depen- 
dent, and is the slope at the faint end of a flux-limited 
survey. Here, k is the lensing convergence, it is related to 
the lensing magnification |^4| by |^4| = 1 + 2k in the weak 
lensing limit. It can be viewed as a projected overdensity 
(e.g. Kaiser & Squires 1993): 



f Xs d\T 

<0,Xg)= —W l (xl,x 9 )5(0,Xl) 
Jo a 



(31) 



where W l ( X l, X 9 ) = (3/2c 2 )H 2 n m r(x L )r(x g - Xl)A(x 3 ) 
is the lensing efficiency function. 

Eq. ( p0| ) tells us how the observed number density of 
galaxies is modulated by gravitational lensing: magnifica- 
tion allows more galaxies to be observed in a flux-limited 
sample, but also causes the galaxies to appear more spread 
out. Which effect wins depends on the slope of the lumi- 
nosity function. 

Gravitational lensing, of course, also modifies the ob- 
served galaxy ellipticity: ei = e- n - + ef x ' (eq. ||). We have 
shown in §[Oj that (8 g e\ n ') vanishes according to two dif- 
ferent tidal theories of intrinsic alignment. Here, we will 
work out (5 g e° x ), for i = t (correlation vanishes for i = r, 
by symmetry under parity). 7 ^| The lensing induced el- 
lipticity for a galaxy at angular position 6 and distance Xg 
away is given by 

ef{0,x g ) = Df [*' -*-.W L {xL,X 9 ri{0,XL) (32) 
Jo a 

where Df = V 2 X - V 2 2 , D e 2 = 2V e iV e 2, and ^(0,XzO = 
20(0, Xi)/r(xi) 2 / (AnGpa 2 ) ((j) is the gravitational poten- 
tial that satisfies the Poisson equation as given in eq. ||To|| ). 

It is important to note that our definition of e° x ' is higher 
than what is usually known as the shear, 7$, by a factor 
of 2. (The usual convention is that if n = \7gip/2 where 
ip is some projected potential, then 71 = (V 2 . — V 2 , 2 ) , 0/2, 
and 72 = V^V^V'-) The reason for our choice is that 
for simple estimators for the quadrupole moment Qij , the 
ellipticity as defined in eq. (|l|) is influenced by lensing 
according to = e" 1 - + e° x - = e\ n - + 27,, at least to the 
lowest order in ellipticity and shear (see Kaiser & Squires 
1993). 

Combining eq. (|), @, (|l|) and @, together with 
the definition of tangential ellipticity in eq. (^|) , we obtain 

(8 g {0)eT-{0')) = (6 g (0)er(6')) A + (6 B {0)<? id')) B (33) 
' J a 2 r{x.hY 

r°° idi 

I ^nk = t/r{xL))J2(t\o-e>\) 



{8 g {d)eT-{d')) B =2 



dXLW b L (xL)W g ( Xg =XL) 



a r(xL) 2 

7 Since, as we will see, (5 9 e^ x ) is non-zero in general, we will 
not consider in this paper the smaller source clustering correction: 
(<5 g e^ ox ) (which is also often ignored in the case of the ellipticity- 
cllipticity correlation). 

8 Strictly speaking, one should consider an additional cross- 
correlation between the magnification bias term involving n and the 
intrinsic ellipticity ej n \ Since n is simply a projected overdensity <5, 
the same arguments in § |2,l| apply directly to this cross-correlation 
as well. 
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f £d£ 

J ■£ p (k = £/r(xL))Mi\o-e'\) 



where 



W[( X l) 



W l (xl) 



dx a W g ( Xg )W L ( XL ,X g ) (34) 



d Xs 5( S - OA)W g ( Xg )W L ( XL , Xg ) 



WKxl) = b / d Xg W g ( Xg )W L ( X L,X g ) 

J XL 

with P{k) being the mass power spectrum, and J% is the 
second order Bessel function: 



J n (y) 



1 

2^ 



d( cos [y sin £ — n£] 



(35) 



Note that W g ( Xg ) is the selection function defined in eq. 
(H) and (30), and Wl{xl) is the lensing efficiency function 
as given in eq. (|3l|). 

We have separated the density-ellipticity correlation 
into two terms: Term A (S g e^')A arises from magnifi- 
cation bias; term B {8 g e c t x ')B arises even if magnifica- 
tion bias is absent. The latter is very much analogous 
to the density-ellipticity correlation commonly measured 
in galaxy-galaxy lensing, except that here we do not mea- 
sure the projected density from one population of galaxies 
and the ellipticity from another - instead, we only have a 
single galaxy distribution W g ( Xg ), from which both den- 
sity and ellipticity is inferred. Note that term B vanishes if 
W g { Xg ) has zero width - it is a non-zero width that allows 
galaxies closer to us to lens galaxies further away from us, 
both sets being drawn from the same W g { Xg ). 

As an illustration, Fig. [l] shows the two contributions 
to (S g et) as a function of angular separation AO (solid and 
dotted lines) for a ACDM cosmological model (il m = 0.3, 
fl A = 0.7, T = 0.21, cr 8 = 0.9), with the following selection 
function: 



W g ( Xg ) = Nz 2 ex V [~(z/0.8y- b ]\dz/d Xg \ 



(36) 



where N is a normalizing factor chosen such that J d Xg W g 
( X g) — !■ Also shown in the figure are the two components 
of the ellipticity-ellipticity lensing correlation (short and 
long-dashed lines): 



(er(0)er(0')> = 4y 



d XL w L { XL y 

a 2 r( XL ) 2 
f°° £d£ 1 

Jo 2lT 1 

d XL W L ( XL ) 2 



{^■{ey*-{e')) 



m 

2vr 



P{k 



t{ X l) 



£/r(x L ))-[J o (iA0) 



(37) 
J 4 (£A9)} 

J 4 (£A0)} 



(Kaiser 1992; Jain & Seljak 1997). The power spectrum 
used is the nonlinear power spectrum appropriate for the 
ACDM model (Hamilton et al. 1991, Jain, Mo & White 
1995, Peacock & Dodds 1996). The model and selec- 
tion function chosen here is in rough accord with exist- 
ing ellipticity-ellipticity measurements from deep lensing 



surveys (e.g. van Waerbeke et al. 2000). The bias pa- 
rameter and the slope of the luminosity function are cho- 
sen to be b — 1 and s = 1, for simplicity. If other 
values were chosen, it is simple to rescale our results: 
(S g e t )A -> (<5 ff e t )ytx(s-0.4)/0.6and (S g e t ) B -» (<$ g e t ) B x& 
(more generally, if b and s were redshift dependent, that 
dependence has to be explicitly integrated; see eq. JUj). 
Note that if s < 0.4, the sign of the the A term is flipped. 
We should also note that on sufficiently small scales e.g. 
AO less than a few arcminutes, the assumption of a linear 
bias likely breaks down; what is shown below nonetheless 
offers a rough estimate of the size of the signal on such 
scales. 
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Fig. 1. — The lensing density-ellipticity correlation (A 
and £?, the magnification bias term and the 'galaxy-galaxy 
lensing' term respectively, see eq. |33[1), and the lensing 
ellipticity-ellipticity correlations (eq. 1371 ) , for a high red- 
shift sample of source galaxies (see eq. [ [36|1 for the selection 
function). 

From Fig. [j], one can see that the A and B contributions 
to the lensing density-ellipticity correlation are compara- 
ble, which are also similar in magnitude to the ellipticity- 
ellipticity correlations. Such a conclusion, however, is sen- 
sitive to the selection function of one's survey. Fig. 
shows the same set of quantities for a source sample with 
lower redshifts: 



W g cx z^exp [-{z/Q3S) 1 - 7 ]\dz/dXg\ 



(38) 



Lower source redshifts decrease the lensing efficiency Wl 
(or Wl', eq. [0), which lead to a drop in all the lensing 
correlations, except the B term for the density-ellipticity 
correlation. This is the term that arises from galaxy- 
galaxy lensing within the sample. Two opposing effects 
roughly cancel out each other in this case: on the one 
hand, lower source redshifts lead to less efficient lensing; 
on the other hand, a shallower survey helps prevent pro- 



jection from washing out the galaxy density fluctuation 
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Fig. 2. — Same as Fig. [I, except that the selection function 
is that of a shallower survey (eq. Pq]). 

The 'galaxy-galaxy lensing' term (B term) can be made 
much smaller if we choose a selection function that is suffi- 
ciently narrow. This is illustrated in Fig. [j] which employs 
the following selection function 



W g cx exp[-(z- if /2/0m 2 ]\dz/dx g \ 



(39) 



Here, the B term is about a factor of 5 smaller than the 
A term, in contrast with the previous two examples. The 
important point to keep in mind is that the magnification 
bias term (A term) and the ellipticity correlation func- 
tions are not sensitive to the width of the selection func- 
tion, because for these quantities, W g enters only through 
convolution with a rather broad lensing efficiency function 
(eq. HH). The B term, on the other hand, can be made 
as small as one wishes if one has the ability to measure 
redshifts accurately. 

3. WHAT ABOUT NON-GAUSSIANITY? 



Are we to conclude from §2.1 that the intrinsic density - 



ellipticity correlation should strictly vanish? A crucial as 
sumption made in the c alcu lation of the intrinsic density 
ellipticity correlation in 



2.1 



is that the fluctuation in grav- 
itational potential ef> is Gaussian random. While this is ap- 
proximately true on large scales, it cannot be exact - grav- 
itational instability induces non-Gaussianity from initially 
Gaussian conditions, even on large scales. It is important 
that we estimate the strength of the non-Gaussian correc- 
tions. To do so properly, we have to make a digression 
and discuss ellipticity-ellipticity correlation. Since non- 
Gaussian fluctuations are easier to include in the frame- 
work of MWK, that is what we are going to adopt. 
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Fig. 3. — Same as Fig. [j] except that the selection function 
is chosen to be much more narrow (eq. [p9[). 



In §3.4, we will also briefly examine the consequence of 
allowing for a nonlinear galaxy- bias. It turns out doing so 
brings in terms quite similar to those introduced by non- 
Gaussianity, as far as the density-cllipticity correlation is 
concerned. Because of the wider relevance of non-Gaussian 
terms (for ellipticity-ellipticity as well as density-ellipticity 
correlations), they will take up most of this section. 

3.1. A Digression on Intrinsic Ellipticity- Ellipticity 
Correlation 

We take as starting point the expression for ef 1 ' given in 
eq. (|jj|), except that we now allow 4> to be non-Gaussian 
in such a way consistent with gravitational instability. 
The usual Gaussian assumption follows from working out 
the angular momentum of galaxies due to tidal torque-up 
according to first order Lagrangian perturbation theory 
(Doroshkevich 1970). A higher order calculation^ would 
obviously introduce non-Gaussian terms, but it should be 
emphasized that a consistent expansion does not necessar- 
ily yield an expression for ef ' that is given exactly by eq. 
(M) where <fi is simply replaced by its non-linear counter- 
part. Nonetheless, this should provide a useful order of 
magnitude estimate. It can be shown that a second order 
expansion does give rise to terms similar to the ones we 
consider here, albeit with different coefficients of the same 
order. Further details will be given in a separate paper 
(Zhang & Hui 2002). 

Because we are interested in the relative significance of 
non-Gaussian versus Gaussian terms, some care must be 
taken in choosing the time at which the spins of galaxies 
(and therefore the tidal fields) are to be evaluated. Ac- 
cording to the usual tidal torque-up theories, galaxies or 

9 See Peebles (1969) and White (1984) for second order calcula- 
tions that assume a spherical Eulerian/Lagrangian volume. 
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halos gain most of their angular momentum before turn- 
around, after which the intrinsic correlations are presum- 
ably locked in place. However, as Porciani et al. (2001a, b) 
pointed out (see also Sugerman et al. 2000), random kicks 
by accreting satellites cause a misalignment of the even- 
tual angular momentum from the one predicted by tidal 
torque-up theories (van den Bosch et al. 2002 also showed 
that hydrodynamic effects cause a misalignment between 
the angular momentum of the gas with respect to that of 
the dark matter). The hope is that a reduction of the 
intrinsic correlations due to such effects occurs in a scale- 
independent way, and can therefore be taken care of by 
rescaling the overall amplitude of the intrinsic correlations 
(as is done by MWK). Porciani et al. (2001a) showed that 
angular momentum growth slows down after z ~ 3, while 
misalignment starts to grow quite early on, from z ~ 50 
onward. In what follows, to illustrate the range of possi- 
bilities, we will therefore consider two cases: the tidal field 
in the expressions for the intrinsic ellipticities (eq. 
are to be evaluated at zt = 3 or zt — 50. p*| 

Let us focus on (^'(Q)^'^')) f° r an order of magni- 
tude estimate: 

<4 n W4 n (0')) = / d Xg W( Xg ) 2 J dp 3 ^ (40) 

([4033012 - 20H012 - 2022012 ~ 60i 3 023]p 
[4033012 - 20H012 - 2022012 ~ 601 3 023]q) 

where we have employed Limber's approximation and used 
the same notation as in eq. (11): pj_ — q± = r(x g )(0 — 8'); 
P is as given in eq. ([)]). 

Let us consider for example the term: 

(033(p)012(p)033(q)012(q)) = (41) 

(033 (p)0i2 (p)) (033 (q)0i2 (q)> 
+ (033 (p)033 (q)) (0i2 (p)0i2 (q)) 

+ (033 (P)012 (q)> (033 (P)012 (q)> 
+ (033(p)012(p)033(q)012(q))c 

The first term on the right can be ignored due to isotropy. 
The next two terms are Gaussian terms which depend 
quadratically on the power spectrum, and are included 
in the analysis of MWK. The last term on the right is the 
non-Gaussian term. One can write the above as: 



{AnGpa 2 
d 3 k A 

d 3 k B 



1 (033 (p)0i2 (p)033 (q)0i2 (q)) 



(42) 



P(fc y 



(2tt) 3 

d 3 k A 



{k A ) A 



■(p-q) 



Bv,B\1 



-ik B (p-q) 



(k B ) 4 



P(k J 



(v.A\2uAi,A 
l.ft-3 ) «■! «-2 



(2tt) 3 ' v " ' (k A Y 
d 3 k A d 3 k B d 3 k c 



(p-q) 



T(k A ,k B ,k c ,k D ) 



(2tt) 3 (2tt) 3 (2tt) 3 " 

(UA\2 h.B h B (h.C\2 hD h D 

(k A ) 2 {k B y (key {k D y 

10 The mapping from Lagrangian to Eulerian space which moves 
the galaxies around also causes an evolution of the intrinsic correla- 
tions. A glance at Fig. 11 of Porciani et al. (2001a) suggests this is 
a small effect. 



where k A + k B + k c + k D = 0, and T is the trispectrum. 
Here, the power spectrum P and the trispectrum refers 
to those of the mass 8, and the factor of (AirGpa 2 )^ 4 (to- 
gether with factors of k 2 ) takes care of the scaling between 
and 8. 

The trispectrum T scales like the power spectrum cubed, 
a scaling motivated by perturbation theory (under Gaus- 
sian initial conditions) and confirmed by numerical N-body 
simulations even in the nonlinear regime (see Fry 1984, 
Scoccimarro et al. 1998, Baugh et al. 1995, but see also 
Suto & Matsubara 1994, Jing & Boerner 1998): 

T(k A ,k B ,k c ,k D ) = (43) 
R a [P(k A )P(k A + k B )P(fc c ) + 11 other cyclic perm.] 
+R b [P(k A )P(k B )P(k c ) + 3 other cyclic perm.] 

The coefficients R a and Rb are configuration dependent on 
large scales, but asymptotes to constants on small scales. 
Both are expected to vary somewhat slowly with scale 
when averaged over configurations (e.g. see variation of 
the kurtosis, which is related to R a and Rb in an averaged 
sense, as a function of scale in Baugh et al. 1995). As a 
rough order of magnitude estimate, we will follow Scoc- 
cimarro et al. (1999) to treat R a ~ Rb ~ Qi, where Q4 
is given by Hyper-Extended Perturbation Theory (HEPT; 
Scoccimarro & Frieman 1999). We will use Q4 ~ 3. One 
should keep in mind that HEPT strictly applies only at 
highly nonlinear scales, and likely leads to some overesti- 
mate of the trispectrum term, but is probably acceptable 
in the spirit of an order of magnitude estimate. 

Putting everything together, one obtains the following 
approximate form for eq. (MX'. 



(44) 



(4 n -(0)4 n -(0')}~[x*y d Xg W( Xg f] 
A ^H^A0 {mmH2o/2r 

10 x* 

(1 + z t ) 4 [ttA 4 (^) + 7rQ 4 A 2 (fc*)A 4 (fe ii ) 
+2^g 4 A 4 (fc,)A 2 (fc i? ) + 7rQ 4 A 6 (fc,)] ZT 

where x* is the typical (comoving) radial distance of the 
galaxies, k* = l/r(x*)/A0, k^ = 1/R where R is the ap- 
propriate smoothing scale for the galaxies (taken to be 1 
Mpc/h by MWK), and A 2 (fc) = 47rfc 3 P(A;)/(27r) 3 , with 
P(k) evaluated at the redshift zt (which will be taken to 
be 3 or 50, higher than the typical redshift of the galax- 
ies which we will refer to as z* (the latter corresponds 
to x*)). We have used the Poisson equation here to 
relate mass fluctuation 8 and (eq. |l0|]), hence the fac- 
tors of 37JQf2 m /2. It should be noted that the integrals 
over multiple power spectra leading to terms like A 4 (fc„) 
or A 2 (fc*)A 4 (fcii) do not strictly factorize, and eq. (44) 
should be viewed as a rough approximation. To be con- 

11 Ours is a particularly demanding application of HEPT, because 
we are evaluating the four-point function at two pairs of points where 
each pair consists of closely separated points, while the two pairs are 
widely separated from each other. 

12 The expression in eq. (|44[), even if one ignores the non-Gaussian 
terms, differs slightly from the corresponding one given by MWK in 
its redshift (2*) dependence. Essentially, MWK sets = z t . To 
the extent that (1 + z) 4 A 4 (fc,z) is only weakly redshift dependent, 
this difference is small. 
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Crete: for instance, one of the Gaussian terms goes like 
d 3 k A d 3 k B 



(2tt) 3 (2tt) 3 



P(k Jl )P(k B )2n6(k A + fcf ) 



G(k A ,k B )e- l ( kA+kB HP-^ 

where the delta function comes from the Limber's ap- 
proximation integral over dpz (eq. p(J), and G repre- 
sents some function of the unit vectors k A = k A /k A and 
k B = k B /fc B . We approximate the above as oc A 4 (A;*)//c*, 
where fc* = l/(p - q)j_ = l/A0/r(x*)- 

For the normalization factor 0, MWK recommended us- 
ing 







1125 0.4 



32 2.8 2 (m m H 2 /2)< 



(45) 



We have introduced an additional overall factor of A 2 in 
eq. (Q), because we need to rescale our result so that 
the correlation still has the desired amplitude after the 
non-Gaussian terms are introduced (also because we have 
chosen to evaluate the tidal field correlations at zt). We 
adopt the choice that (e^'iO)^" 1 '^')) agrees roughly with 
MWK today at fc* = kn i.e. we choose 



A 2 = 



7rA 4 (A: fl ,z = 0) 



(1 + z T ) 4 kA 4 (fc fl ) + 4^Q 4 A 6 (fc fl )] 2=2 



(46) 







(ignoring derivatives). Therefore we 



A crude way to understand our result in eq. (|44|) is this 
Think of ej n - oc 

have (e- n (p)e- n -(q)) very roughly given by 

(0 2 (p)0 2 (q))^(0(p)^(q)) 2 + (47) 

(<Hp)0(q))(0 2 )(0 2 } + 
<0(p)0(q)) 2 (^) + (0(p)0(q)) 3 

where we have not been careful about the coefficients of 
each term. We have also ignored a term that goes like 
(4> 2 ) 2 , which if worked out properly with the right combi- 
nation of derivatives, gives (ef ') 2 = 0. The first term on 
the right is the usual Gaussian result, giving a quadratic 
scaling with the two-point function. The second, third and 
fourth terms arise from non-Gaussianity (using the hier- 
archical scaling of the trispectrum), and the second term 
in particular gives us a linear scaling with the two-point 
function. It should be kept in mind that the above fac- 
torization is actually not exact once derivatives are taken 
into account. One might also worry that the factors of (4> 2 ) 
might yield zero once the various derivatives are properly 
combined. It can be shown that this does not happen. 

The key implication of eq. ( f47j ) (and eq. [Q) is that 
at least one of the non-Gaussian terms is non-negligible 
on sufficiently large scales. At a sufficiently large angular 

13 It shouldie noted, as MWK pointed out, the formulation laid 
down in eq. (bh for the relation between the intrinsic ellipticity and 
the tidal field is, strictly speaking, a little problematic, since it allows 
the possibility of ej n - > 1. They found that in practice the probability 
of that happening remains small as long as they choose f3 properly. 
The tuning of A does the same thing for us here. But one should 
keep in mind a non-Gaussian tail probably makes the problem a little 
worse. 



separation, the non-Gaussian term that scales linearly with 
the two-point function or power spectrum (oc A 2 (fc*)) will 
dominate over the usual Gaussian term (oc A (&»)). More 
precisely, on sufficiently large scales, the dominant term 
scales with the angular separation as A6'A 2 (fc» oc 1/A9). 
At what scale this happens depends on the depth of the 
survey as well as on zt, the redshift at which the tidal 
correlations are to be evaluated. Examples are shown in 
Fig. |j. We use the linear power spec trum corresponding 
to the ACDM model discussed in §2.2 when evaluating eq. 
m4) (and similarly for the rest of this section). 
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Fig. 4. — Ratio of the full expression (including non- 
Gaussian terms; eq. for the intrinsic ellipticity- 
ellipticity correlation to its Gaussian version (i.e. setting 
Q4 — in both eq. |Q and (46)). The upper panel as- 
sumes a selection function appropriate for a deep survey 
(eq. j36|), while the lower panel uses that for a shallow sur- 
vey (eq. |Q). Dotted lines denote the ratio when using 
zt = 3, while solid lines use zt = 50. 

The results in Fig. El might seem a little surprising, es- 
pecially for the case where zt = 50, so that A(fc^) — 0.056 
(recall that A(fc) is evaluated at Zt), which may lead one 
to think that non-Gaussianity can be completely ignored. 
This is true of the normalizing factor (eq. p6[), where 
one can simply ignore the term proportional to Q4. How- 
ever, the same cannot be done for eq. (Q). At sufficiently 
large scales (small fc»), A(fc*) is small enough such that 
A 4 (fc.) < Q 4 A 2 (fc*)A 4 (fc fl ). Q 

14 One might wonder if the significance of the non-Gaussian terms 
can be reduced if the ellipticity is smoothed over some large scale 
by hand. This corresponds to smoothing e^ n - ~ (f> 2 . Because such 
an 'after-the-fact' smoothing is not applied to <f> but to <f> 2 , it can be 
shown that the significance of the non-Gaussian terms compared to 
the Gaussian one cannot be reduced in this way. Therefore, as long 
as (p can be viewed as being smoothed on a small (galactic) scale 
kn ~ 1 h/Mpc, our conclusions remain valid. 

15 We note that in the case of = 3, perturbation theory is close 
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Our argument here is actually quite similar to the well- 
known derivation of linear biasing in the case of the galaxy 
correlation. Suppose S g is some local nonlinear func- 
tion of 5, one can argue on quite general grounds that 
(S g (p)Sg(q)) scales linearly with (J(p)J(q)) for sufficiently 
large |p — q| (Scherrer & Weinberg 1998). 

Since the extrinsic correlation is also expected to scale 
linearly with the power spectrum, our calculation implies 
that the ratio of the intrinsic signal to extrinsic signal does 
not actually drop with scale (on large scales), unlike the 
case of MWK or CNPT who considered the Gaussian term 
only. This is reminiscent of the conclusion drawn by Cate- 
lan et al. (2001), although the reason is quite different - we 
assume here that the galaxy ellipticity depends quadrati- 
cally on the gravitational potential tj> (following MWK and 
CNPT), while Catelan et al. employed a linear relation. 

3.2. Non-Gaussian Intrinsic Density- Ellipticity 
Correlation 

Using very similar arguments to the above, we can write 
down the following approximate expression for the intrin- 
sic density-ellipticity correlation: 

{5 g {ey?-{8')) ~ lx*j dx g W g ( Xg ) 2 }^ (48) 

D(zt) 2 x* 
(1 + z T ) 2 [A 4 (fc») + 2A 2 (fc,)A 2 (fcfl)] 2T 

where we have used the hierarchical scaling: the mass bis- 
pectrumS(k A ,k B ,k c ) = Q 3 [P(k A )P(k B )+P(k B )P(k c )^ 
P(k c )P{k A )] 1 with Q 3 ~ 2 ( using HEPT as before; Scoc- 
cimarro & Frieman 1999). We have ignored the source- 
clustering correction here.[£] The b here is the linear bias 
factor (eq. pi). The factor D(z*)/D(zt) stands for the 
ratio of the linear growth between z* and Zt- This ac- 
counts for the linear growth of S g (or 6) between the two 
different rcdshifts. 

It is interesting to compare this against the extrinsic 
density-ellipticity correlation, e.g. the B-term (eq. Q): 

(8 g (6)(? -)b ~ 3&O m (l + z*)irr( X *) 2 H$A6A 2 (K) (49) 



dXLWg(xL) I Wg( X g)r(xg -XL)/r(x g )} 



where A 2 (fc») is evaluated at z*, the typical redshift of 
the galaxies. The ratio between the intrinsic and extrin- 
sic density-ellipticity correlations is shown in the upper 
panel of Fig. ^|, for a deep survey (z„ ~ 1, W g as given 
by eq. |3(|; long-dashed line for zt = 3, and dotted line 
for zt = 50, where zt is the redshift at which tidal cor- 
relations are evaluated), and a shallow survey (z* ~ 0.3, 
W g given by eq. p8j; short-dashed line for zt = 3, and 
solid line for zt = 50) . Intrinsic alignment induces a non- 
negligible contribution to the observed density-ellipticity 

to breaking down, because A(fc^) approaches unity. 

ie Source-clustering corrections in principle can introduce terms 
that are important as well, but they will not alter greatly our order 
of magnitude estimates, and will not change our conclusion regard- 
ing the scaling of the correlations. We therefore ignore them for 
simplicity. 



correlation, especially for a shallow survey. For a suf- 
ficiently deep survey, however, intrinsic alignment con- 
tributes only at the level of a few percent, if the red- 
shift at which tidal correlations are evaluated, zt, is large. 
The latter suggests that a measurement of both ellipticity- 
cllipticity and density-ellipticity correlations from a deep 
survey might offer a nice consistency check of the lensing 
hypothesis, provided the galaxy bias can be independently 
constrained. We will come back to this point in §[|. 

We caution that our estimates in Fig. on small angular 
scales (less than several arcminutes) should be viewed with 
some skepticism, since 5 (or S g ) has likely gone nonlinear 
in this regime. 
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Fig. 5. — Ratio of intrinsic signal to extrinsic signal. Up- 
per panel: the ratio of eq. ( |48| ) to eq. ( |49| ) for the density- 
ellipticity correlation, for two different selection functions. 
The short-dashed and solid lines correspond to a shallow 
survey (W g given by eq. p8f), with the short-dashed line 
using zt = 3, and the solid line using zt = 50. The 
long-dashed and dotted lines correspond to a deep survey 
(Wg given by eq. Q), with the long-dashed line using 
zt = 3, and the dotted line using zt = 50. Lower panel: 
the level of contamination from intrinsic alignment to a 
galaxy-galaxy lensing measurement of density-ellipticity 
correlation (ratio of eq. to |nj), using two differ- 

ent prescriptions for zt- The dashealine uses zt = 3 and 
the solid line uses zt = 50. The foreground and back- 
ground selection functions here are chosen to mimic the 
SDSS galaxy-galaxy lensing survey of Fischer et al. (2000) . 



3.3. Application to Galaxy-Galaxy Lensing 

The finding that the intrinsic density-ellipticity corre- 
lation is generally non-zero naturally raises the question 
of how this might impact measurements of galaxy-galaxy 
lensing. Typically foreground and background populations 
of galaxies are defined using a simple magnitude cut or 
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photometric redshifts (or in the case of McKay et al. 2001, 
using spectroscopic redshifts to define foreground, but not 
the background) . Generally, the two populations will have 
a non-zero overlap, and the question is whether intrinsic 
alignment constitutes a significant contaminant through 
this overlap. 

The intrinsic density-ellipticity correlation for such a 
set-up is very similar to that given in eq. (fl8|), except that 
one of the W g is replaced by WJ, the foreground selection 
function (the other W g is used to refer to the background) : 

{5 g {eyr{0')) ~ be* J d Xg w g (x g )w g f(x g )} (so) 

— (3Q m H /2) -Q 3 _-_&— 
(1 + z T ) 2 [A 4 (M + 2A 2 (k*)A 2 (k R )} ZT 

Similarly, the extrinsic density-ellipticity correlation (galaxy- 
galaxy lensing) is given by the generalization of eq. (|49|): 

(6 g (0)et x -) ~ 3bn m (l + z«)irr{ X *) 2 Hv 2 AdA 2 (K) (51) 

[ J dxLW^XL) J™ W a (xMx 9 -XL)/r{ Xg )] 

where A 2 (A;*) is evaluated at z*. 

As an illustration, for the background selection function, 
we adopt W g given by eq. (|3|); for the foreground, we use 

WfbCg) = z 2 exp[-(z/0.17) 2 - 3 ]\dz/d Xg \ (52) 

These parameters are supposed to describe the galaxy- 
galaxy lensing configuration of Fischer et al. (2000) (see 
Guzik & Seljak 2001). The ratio between eq. (M) and 
( |5l| ) is shown in the lower panel of Fig. [|: solid fine for 
zt = 50, and dashed line for zt — 3. Therefore, at a scale 
of 10 arcminutes, a conservative estimate of the intrinsic 
alignment contamination to galaxy-galaxy lensing is about 
10%. We will discuss uncertainties in this estimate in §|4|. 

3.4. Nonlinear Bias 

We have assumed a linear galaxy-bias S g cx S so far. On 
small scales, it is almost for certain that the relationship 
between S g and S will be more complicated. Recent work 
using halo models provides perhaps the most sophisticated 
description of the galaxy-mass relation (e.g. Peacock & 
Smith 2000, Seljak 2000, Ma & Fry 2000, Scoccimarro et 
al. 2001). It is in principle possible to use such models 
to predict the intrinsic density-ellipticity correlation. This 
is more properly dealt with in a separate paper. Here we 
would like to simply point out that a nonlinear bias brings 
in terms rather similar to those we have considered. 

Suppose 5 g = bS + b 2 S 2 (see Fry & Gaztanaga 1993), 
where 6 2 is some constant. Then, (5 g (p)e- n -(q)) contains a 
term that goes like b2(S(p)5(p)<f)(q)(j)(q)) , where we have 
suppressed the derivatives on <j). For a Gaussian ran- 
dom field (5 and <fi), this would produce terms that scale 
quadratically with the two-point function. If i5 and <f> were 
non-Gaussian, there can be terms that scale linearly with 
the two-point correlation as well. It is also worth noting 
that on large scales, it is quite likely 62 is small (see e.g. 
Scoccimarro et al. 2001). 



4. DISCUSSION 

Our findings are summarized as follows: 

1. According to two different tidal alignment theo- 
ries (Crittenden et al. 2001 [CNPT] & Mackey et al. 
2001 [MWK]), the intrinsic density-ellipticity correlation 
should vanish exactly. This includes both the lowest order 
term (S g (9)e 1 l n -(9')}, as well as the source clustering cor- 
rection (<5 ff (0)ef in - (9')), where i = 1,2 or i = t (that the 
i = r component vanishes is guaranteed by parity invari- 
ance alone; see eq. Q|). This calculation makes two main 
assumptions: a Gaussian random gravitational potential 
field (as is assumed by CNPT and MWK) and linear bias- 
ing (eq. jn|). While linear biasing is a good approxima- 
tion on large scales, Gaussianity might not be a good one 
even on large scales (see below). 

2. We have computed the extrinsic density-ellipticity 
correlation due to weak gravitational lensing. Only the 
lowest order term (5 g (9)e1 x -(9')) is considered here (i.e. 
the source clustering correction is ignored); it is non-zero 
in general. There are two contributions to it (eq. [j33|) : one 
arises from magnification bias (that gravitational lensing 
modulate the observed density of sources in a flux limited 
survey; this we call the A-term), and the other is analo- 
gous to what is commonly known as galaxy-galaxy lens- 
ing, except that here we are interested in galaxies drawn 
from a single selection function (foreground galaxies lens- 
ing background galaxies from the same selection function 
of some finite width; this we call the i?-term). Illustra- 
tions of these two terms, compared against the more fa- 
miliar ellipticity-ellipticity lensing correlations, are shown 
in Fig. [l]to |[ Two main trends are noteworthy. First, the 
B term is systematically higher than all other terms if the 
source galaxies are at sufficiently low redshifts (z < 0.3); 
by z ~ 1, all terms are roughly comparable. Second, the A 
term can be made dominant over the B term if the width 
of the selection function can be made sufficiently narrow. 

3. We have considered non-Gaussian contributions to 
both ellipticity-ellipticity and density-ellipticity intrinsic 
correlations. They are non-negligible. 

a. In the case of the intrinsic ellipticity-ellipticity cor- 
relation, the non-Gaussian contributions predict a large 
scale linear scaling with the power spectrum rather than a 
quadratic one as discussed by MWK and CNPT. The re- 
sult is somewhat surprising, in that higher ord er t erms 
actually dominate over lower order ones (see § |3.l| ). It 
originates from the fact that the lowest order Gaussian 
term scales with the angular separation like A#A 4 (fcA 
while one of the higher order non-Gaussian terms scales 
as A#A 2 (/c*), where A(fc«) is the fluctuation amplitude 
at wavenumber k* oc 1/A9 (eq. Even though the 
higher order term is suppressed by a small coefficient, it 
still dominates at sufficiently large angles A9. The am- 
plitude of the resulting ellipticity-ellipticity correlation is 
somewhat uncertain - we fix it by matching the results of 
MWK on small scales (so as to obtain the right rms value 
for the observed ellipticities) ; but the large scale scaling 
we believe to be robust. 

b. We compute the intrinsic density-ellipticity correla- 
tion, which receive contributions from non-Gaussian terms 
alone. A comparison of this with the extrinsic density- 
ellipticity correlation is shown in Fig. ||. For a shallow 
survey (e.g. median redshift of 0.3), the intrinsic signal is 
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non-negligible. However, for a sufficiently deep survey (i.e. 
median redshift of 1 or above), the intrinsic signal is only 
a small fraction (a few percent) of the extrinsic signal on 
large scales, if Zt, the redshift at which tidal correlations 
are evaluated, is large enough. 

c. We apply the calculation of density-ellipticity correla- 
tion to the case of galaxy-galaxy lensing, and find that for 
a low redshift survey such as SDSS (Fischer et al. 2000), 
intrinsic alignment constitutes a non-negligible contami- 
nant, roughly at the 10 — 30% level at 10 arcminutes (lower 
panel of Fig. ||). 

There are several issues that are worth exploring in the 
future. 

First, we advocate the measurement of density-ellipticity 
correlation from current weak-lensing surveys. This can 
be straightforwardly implemented using current data and 
increases their scientific return. Most surveys are suffi- 
ciently deep that the intrinsic contribution should be small 
( a fe w percent; as long as zt is sufficiently large - see 
§ |3.2| ). The same conclusion was reached regarding the 
eHipticity-ellipticity correlation by MWK and CNPT. It 
would be useful to check that this is indeed the case, by 
measuring both correlations, and see if they are consistent 
with each other under the lensing interpretation. Such a 
check requires knowledge of the galaxy bias b as well as 
the luminosity function (slope s), however. The luminos- 
ity function is in principle directly measurable, while the 
galaxy bias b can be obtained from higher order clustering 
measurements (see Fry 1994, Scoccimarro et al. 2001 and 
references therein) . What makes this program challenging 
is that both s and b can be redshift-dependent. 

Second, it is important to check for the possible contam- 
ination of low redshift galaxy-galaxy lensing measurement 
by intrinsic alignment. Better measurements of intrinsic 
alignment in surveys of closeby galaxies will be quite in- 
teresting (e.g. Brown et al. 2000). The SDSS survey is in 
principle very useful for this purpose (McKay et al. 2002, 
unpublished preprint) . The different scalings with redshift 
between the intrinsic and extrinsic signals can also be used 
to tell them apart (see eq. [|o| & Q). 

It would also be important to check the predictions 
we make in this paper, for both ellipticity-ellipticity and 
density-ellipticity correlations, against numerical simula- 
tions, especially on large scales (above several Mpc/h). 
A non-zero measurement of the density-ellipticity correla- 
tion on large scales from simulations might be interpreted 
as indicating that non-Gaussian tidal fluctuations are in- 
deed important, since Gaussian theories predict a vanish- 
ing contribution. We should emphasize that the normal- 
ization of the correlations is really not predicted by current 
analytic theories. It involves free parameters, such as the 
moments of inertia, which have to be fixed by matching 
simulations and/or observations. Presumably, such a nor- 
malization procedure can approximately account for the 
reduction of alignment correlations seen in simulations, 
due to nonlinear effects (e.g. Porciani et al. 2001a, b), 
or hydrodynamic effects (e.g. van den Bosch et al. 2002). 
But it is also quite possible that a more radical change in 
alignment theories is required; testing this will be impor- 
tant for progress. It is for this reason that our estimates, 
such as the level of contamination to galaxy-galaxy lens- 
ing due to intrinsic alignment, should be viewed with some 
skepticism. If the basic picture adopted in this paper holds 



up, the prediction for a large scale linear scaling of the in- 
trinsic correlations with the mass power spectrum should 
be fairly robust. This has important implications, in that 
the ratio of the large scale intrinsic to extrinsic signals does 
not drop significantly with scale. 

We thank Ue-Li Pen for useful comments and Roman 
Scoccimarro for helpful discussions. Support for this work 
is provided by an Outstanding Junior Investigator Award 
from the DOE and an AST-0098437 grant from the NSF. 

APPENDIX: 

ESTIMATORS FOR THE DENSITY-ELLIPTICITY 
CORRELATION IN REAL AND FOURIER SPACE 

An estimator for the real-space density-ellipticity corre- 
lation is 



Across (A0) 



£ a ,g("(«)-")et(/?)n(/3)e^ 



(53) 



where one can imagine pixclizing the survey so that each 
pixel contains either one or no galaxy, symbolized by n(a) 
which equals 1 or depending on whether the pixel a con- 
tains a galaxy or not; n is the average number of galaxies 
per pixel, and 9^' is equal to 1 if pixel a and pixel (3 are 
separated by AO, zero otherwise. This estimator on the 
average provides a measure of Across in eq. (0). The pix- 
elization discussed above can be thought of as a conceptual 
device. In practice, one way to carry out the estimation 
is: take all pairs of galaxies at a given separation of in- 
terest AO, compute the average e t , and then do the same 
for a pair consisting of a galaxy and a random point; sub- 
tracting the two and multiplying the result by appropriate 
constants will yield Across ■ 

Another method to estimate Across (AO) makes explicit 
use of pixelization. Suppose the survey is pixelized so that 
each pixel contains at least several galaxies. Let N(a) 
be the number of galaxies in pixel a. We define a pixel 
ellipticity by 

N(f3) 

6*09) = J2 *S)/N (54) 

2 = 1 

where N is the average number of galaxies per pixel, and 
i labels individual galaxies within the pixel f3. Note that 
et at pixel (3 is defined with respect to another pixel a (in 
order to define the tangential part). 
Then, £ croS s can be estimated by: 



Across (AO) — 



E a AN(a)-N)e t (f3)Q^ 



(55) 



This estimator would also yield Across (eq. (?]]) on the av- 
erage. 

The reader might wonder why we had not used a seem- 
ingly more natural definition of pixel ellipticity: et(/3) = 

J^^JP e t(i)/N(l3). A moment's thought would reveal that 
using this instead in eq. ( ^5| ) gives 

, / dX g W g {x g ) J dx' g W g ( X ' g )S g (8, Xg )(l+S g (8' , X ' g ))e t (6' , X ' g ) 
{ J dx^W g (x^)(l+S g (9', X ' g ')) ' 



rather than eq. (|7|) . The above is a legitimate quantity to 
consider, but the expectation value then no longer simply 
consists of two terms as in eq. (Q). Rather, there will be 
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an infinite number of terms (see Bernardeau 1998 and Hui 
& Gaztanaga 1999 for discussions of related issues). 

In this paper, we have focused exclusively on the density- 
cllipticity correlation in real space. Let us provide here the 
corresponding estimator in Fourier space. A Fourier space 



estimator has the usual advantage of uncorrelated band- 



p ower errors, the complication ot window function aside. 
We will continue to think of a pixelized survey where each 
pixel has at least several galaxies. 

Let us define 6(a) — (N(a) — N)/N and its Fourier 
transform: 



6(1) =d 2 8j2kc 



t£-0 r 



(56) 



where d 2 9 is the angular size of each pixel. 
Let us also define the Fourier transform: 



~e j (£) = d 2 6T / eMe 



-ite 







/3 



(57) 



where j = 1,2, and £j(j3) = ^2^JP £j(i)/N . Define also 
e t (£) = [(i\ - e 2 )/i 2 ]e x (t) + [2£i£ 2 /e 2 ]e 2 (£) (note that e t is 
not the Fourier transform of et). 

Finally, we can form the following estimator: 



S) 



i Efmw) 



V 



(58) 



where V is the area of the survey, and J^f refers to sum- 
mation over all modes with \£\ = i. Since the intrinsic 
density-ellipticity correlation vanishes, this would on the 
average gives us an estimate of something related to the 
extrinsic density-ellipticity correlation in Fourier space. It 
can be shown that, ignoring source clustering corrections, 

(PcrossW) = 2 / % ^#^ P(^ = e/r( XL )) (59) 



P(k = £/r( XL )) 



a 2 t(xl) 2 

d X L W b L ( XL )W g ( Xg = Xl) 
a r( XL ) 2 



The two terms on the right hand side are Fourier analogs 
of the A and B terms in eq. (p3[). 
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